Theory of nonlinear optical response of ensembles of double quantum dots 
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We study theoretically the time-resolved four-wave mixing (FWM) response of an ensemble of 
pairs of quantum dots undergoing radiative recombination. At short (picosecond) delay times, the 
response signal shows beats that may be dominated by the subensemble of resonant pairs, which 
gives access to the information of the inter-dot coupling. At longer delay times, the decay of the 
FWM signal is governed by two rates which result from the collective interaction between the 
two dots and the radiation modes. The two rates correspond to the subradiant and superradiant 
components in the radiative decay. Coupling between the dots enhances the collective effects and 
makes them observable even when the average energy mismatch between the dots is relatively large. 

PACS numbers: 78.67.Hc,78.47.Cd,03.65.Yz 



I. INTRODUCTION 

Double quantum dots (DQDs) are pairs of quantum 
dots (QDs) placed at a small distance from each other 
and coupled either by long-range dipole forces^ or by 
carrier tunneling through the inter-dot barrier^^ (in the 
latter case, the system is referred to as a quantum dot 
molecule). Such systems attract much attention both 
in the theoretical and experimental research. This still 
growing interest is driven to a large extent by the techno- 
logical promise these structures show for nanoelectronics 
and quantum information processing applications. A sys- 
tem built of two QDs can be viewed as a first step towards 
scalable semiconductor-based quantum devices. Cou- 
pling between the dots is essential here, as it allows one 
to induce conditional dynamics in the system and thus to 
realize the basic elements of quantum computing 6 . An- 
other interesting application of coupled QDs is generation 
of entangled photons 7 . It is therefore not surprising that 
much experimental work has been devoted to proving the 
existence of coupling in DQD systems ^ 10 ' 11 ' 12 . 

In this paper, we study the four-wave mixing (FWM) 
optical response of an inhomogeneous ensemble of self- 
assembled double quantum dots. The systems under con- 
sideration are QD pairs formed spontaneously by strain- 
induced nucleation, as is typical, e.g., in the InAs/GaAs 
syste m 13 ' 14 (although, in fact, the exact arrangement of 
the QDs is not essential in our discussion). The FWM 
spectroscopy is an optical technique commonly used for 
extracting information on lifetimes and homogeneous de- 
phasing from inhomogeneous ensemble a 15 ' 16 ' 17 . It has 
been also applied to DQD ensembles 1 ^, showing features 
that are clearly distinct from those observed in ensem- 
bles of individual QDs, like two-component decay and 
modified initial dephasing. As the physical properties 
of DQDs are much more complex than those of individ- 
ual QDs it is not always possible to identify the mech- 
anism responsible for these differences. In particular, it 
is not clear to which extent they may result from sim- 
ple optical interference or quantum-optical mechanisms 
of collective radiative decay (sub- and superradiance), as 
opposed, e.g., to DQD-specific phonon-related dephasing 



channels, like dissipative exciton transfe r 19 ' 20 . 

Superradiance has been thoroughly studied for 
strongly excited atomic samples 21,22 , where it is mani- 
fested as an outburst of radiation due to constructive in- 
terference of quantum-mechanical amplitudes for radia- 
tive transition in a massively correlated state of atoms in- 
teracting with a common radiation reservoir. Somewhat 
less spectacular form of this effect was also observed in 
QD ensembles^ 3 - where the radiative recombination was 
shown to increase for a sufficiently large number of dots. 
However, the essential features of the collective emission 
can be observed already in a two-emitter system: states 
with a delocalized excitation decay slower or faster, de- 
pending on the relative phase between the two compo- 
nents in the superposition (the first or the second emitter 
excited)^^. 

The goal of the present work is to identify the purely 
optical effects that may appear in the nonlinear optical 
response of an ensemble of DQDs. It has been pointed 
out that the dynamics of double dots coupled to a ra- 
diative reservoir is very reach, both in the open space 2 ^ 
and in a cavity 26 -. Here, we will show that the exponen- 
tial decay of optical coherence, characteristic of a sin- 
gle QD, is replaced by a non-exponential, two-rate decay 
which may be related to sub- and superradiant compo- 
nents in the system evolution. As in the previously stud- 
ied case of a single pair of nearly identical QDs^, the de- 
cay is strongly affected by the coupling between the dots. 
We show that the collective features in the radiative de- 
cay persist even for relatively large values of the energy 
mismatch between the dots (many orders of magnitude 
larger that the emission line width) as long as the cou- 
pling is of comparable magnitude. Another feature that 
emerges from our model is a damped oscillatory behavior 
of the FWM signal at short delays (sub-picosecond and 
picosecond time scales), which results from the optical in- 
terference of the signals emitted by various DQDs in the 
ensemble, dephased due to inhomogeneity of DQD pa- 
rameters. As we will show, these beats may be dominated 
by a relatively long-living contribution from the minor- 
ity subensemble of resonant DQDs (formed by dots with 
identical transition energies). In this way, the nonlinear 
response contains information on the coupling between 
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the dots, irrespective of the energy mismatch between 
them. 

The paper is organized as follows: In sec. Owe present 
the system under consideration and define its model. 
Next, in sec. IIIII we discuss the evolution of the system 
under two-pulse optical driving and derive the equations 
for the third order optical polarization. The results are 
discussed in sec. |IV] Sec. |V] contains concluding re- 
marks. 



II. THE SYSTEM 

We study an ensemble of DQDs, each consisting of two 
QDs. The QDs in each pair are placed at a distance much 
smaller than the relevant photon wavelength so that spa- 
tial dependence of the electromagnetic (EM) field may 
be neglected (the Dicke limit). We assume that the po- 
larization of the laser beam is chosen such that only one 
of the fundamental excitonic transitions in each dot is 
allowed. Each DQD can be then modelled as a four- 
level system with the basis states |0) = |00), |1) = |01), 
2) = 1 10) , |3) = |11), corresponding to the ground state 
(empty dots), an exciton in the second and first QD, and 
excitons in both QDs, respectively. A single DQD is com- 
posed of two QDs with energies = E ± A. We will 
describe its evolution in a frame rotating with the fre- 
quency E/h. Then, in the rotating wave approximation, 
the Hamiltonian for a single DQD is 

H = Hx + Hl + Hrad + -Hphot- 

The first component describes the excitons 
Hx = 

A (n x -n 2 ) + V (V'Vf + <7<V 2) ) + V B nin 2 [l) 



where <7+ are creation and annihilation operators of an 



where 



exciton in the j-th QD, fij = crV 'a_ is the occupation 
of the dot j, V is the coupling between the dots (e.g., 
tunneling or Forster), and Vb is a biexciton shift due to 
a static dipole interaction 2 -. 

The second term in the Hamiltonian accounts for the 
interaction with the laser field, which is treated classi- 
cally, 
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where 4>i, and U are the amplitude envelopes, phases, 
and arrival times of the laser pulses, respectively. 

The third term accounts for the interaction with the 
quantized EM field (radiation reservoir) in the dipole ap- 
proximation. 



-Hrad = (T 



S 1 ) 



(2) 



g k \ =id- e x (k) 



2t e r v ' 



k is a photon wave vector, Uk is the corresponding fre- 
quency, A denotes polarizations, bk,\,b\. x are photon 
creation and annihilation operators, d is the interband 
dipole moment (for simplicity equal for all QDs), e\ (fc) 
is a unit polarization vector, e is the vacuum permit- 
tivity, e r is the dielectric constant of the semiconductor, 
and v is the normalization volume for the EM modes. 
Finally, 



H 



phot 



(3) 



fc.A 



is the Hamiltonian of the radiation reservoir. 

To describe the ensemble of DQDs, we assume a Gaus- 
sian distribution function for the energies of the two dots, 



g(E 1 ,E 2 ) 



1 



27TCT 2 v /r 3 p f 



(Ej -El) 2 -2p(E 1 -E 1 )(E 2 -E 2 ) + (E 2 -E 2 ) 2 



with the mean transition energies E%, E 2l identical en- 
ergy variances a 2 for both QDs, and a correlation coef- 
ficient p. Note that this distribution corresponds to an 
uncorrelated Gaussian distribution of the parameters E 
and A, g (E, A) = g E (E) g A (A), where 



9a(A) 
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2tt<7a 



exp 



2a\ 



A = E,A, (5) 



fc.A 



with the mean values E — (Ei + E 2 )/2 and A = (Ei — 
E 2 )/2 and variances a\ = er 2 (l + p)/2 and a\ = cr 2 (l — 
p)/2 (correlation between the QD energies E\, E 2 means 
less variance of their difference). 



III. THE SYSTEM EVOLUTION AND FWM 
RESPONSE 

A FWM experiment, which we want to model, consists 
in exciting an ensemble of DQDs with two ultrashort laser 
pulses, arriving at t\ = —t and t 2 — 0. The first step 
of the calculation is to find the optical polarization of a 
single DQD after the second pulse, which is proportional 
to 

P(t) = Pw (t) + p 2Q (t) + p 31 {t) + p 32 (t) + c.c, 

where t > 0, p(t) is the density matrix of a DQD struc- 
ture, and pki — (k\p\l). The first two terms are exci- 
ton coherences (polarizations) while the other two are 
referred to as biexciton polarizations. In order to ex- 
tract the FWM polarization we pick out only the terms 
containing the phase factor e^ 2 ^ 2- ^ 1 ', which mimics the 
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frequency shifting and lock-in detection technique used in 
the actual experimen t 15 ! 27 . In the second step, the total 
optical response from the sample is obtained by sum- 
ming up the contributions from individual DQDs with 
the weight factor g(E, A), 

P FWU (t) = J dAdEg(E, A)P(t). 



A. Single DQD evolution 

The detection of weak signals originating from the 
DQD ensemble is based on a heterodyne technique 15 : 
The response Pfwm is superposed onto a reference pulse 

E Ic{ (t) = / ref (t - t )e- m ^/ n + c.c, 
arriving at a time to. We assume a Gaussian envelope 



/ref(i) - 



1 



■ exp 



1 - 

2 Vr ref 



V27TT re f 

The measured signal is proportional to \F(t 0l r)|, where 



F(t Q ,r) = e ^(*o-)A / d tP^ M (t)E^ (t), (6) 



where -fp^M an< ^ ^Vef are * ne positive frequency part of 
the FWM signal and the negative frequency part of the 
reference pulse, respectively, and the (irrelevant) phase 
factor has been extracted for convenience. 

We assume that the pulses are spectrally very broad to 
assure resonance with all the QDs in the ensemble. If the 
durations of the pulses are much shorter than both h/A 
and h/V, the action of each of them corresponds to an 
instantaneous, independent rotation of the state of each 
QD, that is, to the unitary transformation U; = Ui <g> E/j, 
where 



Ui 



(?) 



— I sin 



o, 



-i(<t>i+EU/h)\ 



Here I denotes the identity operator and 



1 

ai = h 



fi(t)dt 



is the pulse area. 

We assume that the initial state of a DQD is p(— t~) = 
|00)(00| denotes just after or before a time instant t). 
The DQD is then excited with the first pulse. Just after 
this pulse, the system state is 



P (-T+) = Ulp(-T")Uj 



(J) 



and the four matrix elements related to optical polariza- 
tions have the values 

Poi(-r-) = p 02 (-r-) = ^sina 1 cos 2 ^e-* +i;B ^, 



Pl3{-T ) = p23(-T ) 



l - sin ai sin 2 ^e"*^/*. 



From now on, we will only keep the terms which are of 
the first order in the first pulse area a\. In this approx- 
imation, the biexciton coherences do not contribute at 
this stage of the evolution. 

Between the laser pulses, the evolution of the reduced 
density matrix of the charge subsystem is described by 
the Lindblad equation of the form 2 ^ 



p = --[H x ,p]+£[p], 



with 



c[p] = r 



(8) 



(9) 



where T is the spontaneous decay rate for an individual 
dot and E_ = = a ( } ] + a {T> . This yields a closed 
system of four equations of motion for the negative fre- 
quency parts of optical polarizations, 

An - {iA/h-T/2)p 01 + (iV/H-T/2)p 02 

+r(pi 3 + p 23 ), (10a) 
P02 = (iV/h-T/2)poi-(iA/K + r/2)p 02 

+r(pi 3 + p 23 ), (10b) 
P23 = {iV B /h+iA/h-3T/2)p 2 3 

-{iV/H + T/2)p 13 , (10c) 
Pis = -(iV/H + T/2)p 23 

+ {iV B /ti-iA/h-3T/2)p 13 . (lOd) 

The solution to these equation simplifies if one notes 
that for typical double dots the energy mismatch is much 
larger than the fundamental line width. Therefore, in the 
following we assume T <C O/fi, where fl = (A 2 + V 2 ) 1 / 2 
corresponds to the half-splitting of the single exciton 
states. This reflects the actual experimental situation 
as the energy mismatch observed in real samples ranges 
from a few meV to tens of me V 18 ' 28 i 29 while the recom- 
bination rate 1/r is typically of the order of 1 n a 18 i 29 ' 30 . 

Upon solving the equations (|10al d). one finds the ex- 
citon coherences at t — 0~, 

Poi(02)(0-) = \ sin ai cos 2 a e -*i-H*T/ft 



1 - 



V ± A 



with 



n i 

A (i) = Tij- ~ ^T r > 



(12) 



where 0± = 1 ± V/O, and the upper sign corresponds to 
the first index or pair of indices. 

The second pulse, with an area a 2 arrives at t = 0. In 
order to find the FWM response to the leading (third) 
order, we need to calculate the positive frequency parts 
of optical polarizations, keeping only terms of the sec- 
ond order, containing the phase factor e 2 ^ 2 . Such terms 
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depend only on the values of the negative frequency po- 
larizations before the pulse and have the form 



Pio(0+) = ,020(0+) 



P31(32)(0 + ) 



; a 2 [poi(0-)+po2(0-)](13a) 



; "2ft)2(01)(0 )• 



(13b) 



Note that biexciton polarizations appear in the leading 
order and cannot be eliminated since, contrary to the 
single dot case, the transition to the molecular biexciton 
(excitons with the same polarization confined in different 
dots) is not forbidden by selection rules for any polariza- 
tion of the laser beam. 

The single exciton polarizations at an arbitrary time 
t > are found by using Eqs. (|13al b) as initial values 
for the system of equations of motion for positive fre- 
quency polarizations, which is obtained from Eqs. (jlOat -d) 
by complex conjugation. The total single-exciton polar- 
ization is PW(t) = pw(t) + p2o(t) and is explicitly given 
by 

P«(i) = -sinaxCOS^sinW^-^e^ 
w g 2 



■j3+e 



Although, in principle, the biexciton term affects the ex- 
citon coherences in the same order of the optical response 
(due to radiative recombination) one finds that the cor- 
responding terms are of the order of fir/fi <C 1 and can 
be neglected. 

For the biexciton polarizations, the evolution equations 
Paid) yield 



P31(32) 
1 

2 
1 

+ 2 



it 



V 



P31(32)(0 + ) + ^P32(31)(0 + ) 



1 ± 



031(32) (0" 1 



P32(31)(0 + ) 



where 



^2(3) 



Tin/h + iv B /h- [i + -/3± r 



(15) 



The biexciton contribution to the coherent polarization 
is then 



P^{t) 



03 1 + P32 

1 sinai cos 2 21 S in 2 aae <(2*i-*0 e ^r/ft 
8 2 



x 



(ft 



A 2 £+Ai r 



(16) 



B. Ensemble response 



The FWM polarization is obtained upon returning to 
the Schrodinger picture, which amounts to inserting the 



phase factor e ~ tEt / h and adding up the contributions 
from all the dots in the ensemble according to their sta- 
tistical distribution given by Eqs. and (|5|). For the 
sake of the further discussion it is convenient to split the 
exciton contribution (|14p into two parts and to write the 
total polarization as a sum of three contributions 



Pfwmw — 

^ sin (ai) cos 2 (^-J sin 2 (a 2 ) 



3 

/ dAdEg(E, A)e-' lE< - t - T ^ h J2 p n(t, A) + c.c 

i . , s 2 f a i\ ■ 2 / \ 
-sm(ai)cos (^j sin (a 2 ) 



xe 



--iE(t-T)/H-^(t + T) 

3 



J dAg(A)Y,P n (t,A)+c.c, 



(17) 



where we inserted the definitions (TT2|) and (fT5|) into 
Eqs. (O and dTBJ) and defined 



P n (f,A) = P n+ (t,A) + P n _(t,A), 



P 1± (i,A) = /%exp 



VT 

T^(r + t)±m (r-t)/h 



P 2± (i,A)=/3 ± £_ exp 
P 3 ±(t,A) = 



VT 

±X}(T-t)T*n (r + t)/h 



exp 



{T + iV B /h)t± 



VT 

m 



in/h) (r + t) 



The second form of Eq. (TIT)) is obtained by performing 
the integration over E. 

Next, one has to calculate the heterodyne signal gener- 
ated by the FWM polarization overlapped with the ref- 
erence signal, according to Eq. ([6]) . The integration over 
time in Eq. ^ can easily be performed. Substituting 
Eq. (TIT)) into Eq. ^ and performing some reasonable 
approximations, as discussed in detail in the Appendix, 
one finds 

3 



where 

F n (t ,T) 

Here 



■ sin (ai) cos 2 sm2 ( a 2) 

J dA 5A (A)[jF„ + (A)+JP n _(A)].(18) 



:F„±(A) - e iE{to-T)/H dt pi+\t,A)E(J(t) 



exp 



°i(*o-r) 



2 1 



2 (ft 2 



'rcf 



x exp 



nfo - r) 



h{h + Tl ( al/h) 



$n±(A), (19) 



5 



where 
*i±(A) 

$ 2± (A) 

$ 3 ±(A) 



/3 ± exp 



2^2 



x exp [— /3±IY] , 



/?+/?- exp 



' rof 



2(ft 2 



x exp [=p2i0r/fi] exp [— rr] , 



(20a) 
(20b) 



exp 



2(h 2 



cf u E) 



x exp 



iV B (t 



exp[-(l + /3 T )rr] 



(20c) 



In order to calculate the time-resolved nonlinear re- 
sponse of the DQD ensemble, Eq. (fT9|) must be inte- 
grated numerically with the distribution function g&(A), 
according to Eq. (fTS|) . Another integration, over to, yields 
the time- integrated (TI) FWM signal as a function of the 
delay time r, which is commonly used to characterize the 
dephasing in a physical system. The results will be dis- 
cussed in the following section. 



IV. RESULTS 

In this section we present and discuss the time-resolved 
and time-integrated FWM response from an ensemble of 
double quantum dots, depending on the statistical distri- 
bution of the energy mismatch between the dots in the 
ensemble and on the strength of the coupling between 
the dots. We assume fixed values of the average energy 
mismatch A = 4 meV, the standard deviation of the 
mean DQD transition energy erg = 8 meV, the length of 
the reference pulse r re f = 21 fs (corresponding to 50 fs 
full width at half maximum), and the spontaneous re- 
combination rate for an individual dot T = 1 ns _1 . We 
start by studying the nonlinear response for short (pi- 
cosecond) delays; then we proceed to the discussion of 
the signal decay on long (nanosecond) time scales. 



A. FWM response for short delay times 



It follows from Eq. ([19]) that all the contributions to 
the FWM signal are restricted to the short range of delay 
times around to — T , of width ~ Ti/cte, that is, they 
have the form of a photon echo. However, among the 
three contributions F n [Eq. (fT8|)]. the first one, F\, has a 
different character than the other two. This results from 
the different structure of the phase factors of the form 
gin-r/fi appearing in the last exponential term of Eq. (flT))) . 
in the second exponential term of Eq. (|20b|) and in the 
second exponential term of Eq. (|20c|) . Since fl depends 
on A, which varies across the ensemble, such terms tend 
to interfere destructively when the signal from different 
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FIG. 1: (Color online) (a-c) Real and imaginary parts (green 
dashed and blue dotted lines, respectively) of the three contri- 
butions F n (t ,r) to the FWM echo F(to,r) [Eqs. lfl8 ]l -pDc ) ] 
for short delays, with r = 0.23 ps, V = 2 meV, Vb = 1 meV, 
<ta — 1 meV. (e-g) As previously but at r = 0.61 ps. (d,h) 
Real and imaginary parts of the total signal F(to,r), as well 
as its modulus, which corresponds to the measured signal, for 
the two values of r. 



DQDs is added up. However, in the case of Ft, this 
phase term depends only on to — t, which is limited to 
the width H/ct-e of the echo pulse. Therefore, the spread 
of the phase factors is also limited and independent of 
t. The only effect of this phase term is therefore a slight 
asymmetry of the echo pulse due to oscillations in Im Ft , 
as shown in Fig. [lja,e) (which makes it different from 
the simple pulse shape from an ensemble of individual 
dots). These oscillations are always in phase with the 
center of the echo peak, so that the peak area (that is, 
the time- integrated signal) is constant. 

The situation is different in the case of the other two 
contributions F% and F3. Here, another phase term ap- 
pears, proportional to f2r jfi or (Vb ± 2Sl)r/H [in the sec- 
ond exponential terms of Eqs. (|20b[) and (|20c|) ]. Thus, 
the phase of the signal at to — r varies with r, which 
leads to a variation of the shape and magnitude of the 
echo. As can be seen by comparing Figs. HJb,c)] with 
Figs. njf,g)] , the contribution of these terms to the pho- 
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1 2 2 4 6 8 

T [ps] t [ps] 

FIG. 2: (Color online) (a,b) The evolution of the time- 
integrated signal for short delays without (a) and with (b) 
a biexciton shift for <ta <C A (no resonant dots in the en- 
semble). Red solid lines: V — 2 meV, <ta = 1 meV; green 
dashed lines: V = 2 meV, <ta = 0.5 meV; blue dotted lines: 
V — meV, tj a = 1 meV. (c,d) As in (a) and (b), respectively, 
but for o a comparable with A. Red solid lines: V — 2 meV, 
<7a = 3 meV; green dashed line: V = 2 meV, <ta = 5 meV; 
blue dotted lines: V = 3 meV, cta = 3 meV; grey dash dotted 
line V = 0, <ta = 3 meV. 



ton echo strongly depends on the delay r. As a result, 
also the amplitude of the total (actually measured) sig- 
nal varies with time on picosecond time scales, as shown 
in Figs. QJd,h). This variation leads to oscillations in the 
TI signal (Fig. [J). These oscillations are a manifesta- 
tion of optical beats between the two dots. Their form 
depends on whether the ensemble contains a fraction of 
DQDs composed of identical (resonant) dots, that is, on 
the interplay of a a and A. 

If a a "C A then the signal originates from all the 
DQDs, whose values of the energy mismatch A lie 
roughly within a a from A. The inhomogeneity of the 
values of A translates into inhomogeneity of £1 and leads 
to a spread of phases in the terms like exp(«fir)/S, which 
increases as r increases. Since the total signal from the 
sample is a coherent sum of the fields emitted by all 
DQDs, this phase distribution leads to quenching of the 
two contributions Fi and F3 at r ~ Ml /(Act a) « 1/r 
and, therefore, to vanishing of the oscillations. This can 
be clearly seen in Fig. [2la,c). This effect is due to the fact 
that the probe pulse acts symmetrically on both dots and, 
therefore, can invert (refocus) only the dephasing due to 
the inhomogeneous distribution of the average transition 
energies but not that resulting from the inhomogeneity 
of the energy mismatch between the dots in a single pair. 

The evolution of the optical signal becomes more inter- 
esting if a a ^ A. This condition means that the ensem- 
ble contains a fraction of resonant DQDs, that is, such 
that have nearly identical transition energies. The fre- 
quency fl has a minimum at A = which corresponds 
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FIG. 3: (Color online) Time-integrated FWM response, (a) 
Dependence on the coupling strength V for cta = 3 meV. (b) 
Dependence on the inhomogeneity of the energy mismatch for 
V = 5 meV. 

to a stationary point of the phase distribution over the 
ensemble. Therefore, all the nearly resonant QDs emit 
radiation in phase and can dominate the contributions F2 
and F3 of the ensemble when the signal from the possibly 
much more numerous dots with A ~ A has dephased as 
described above. An essential point is that, since this mi- 
nority resonant subensemble has A ~ 0, the frequency of 
the beats is very close to 2V/H or (2V — Vb) /fi [note that 
the term with this frequency has a much larger amplitude 
in Eq. (|20cf than that with the frequency (2V + Vb)/H 
since for A — > one has (3- — > 0]. In this way, the non- 
linear response gives a direct access to the values of the 
inter-dot couplings and to the properties of resonant dots, 
even if they are a minority in the ensemble. The beats 
originating from the resonant subensemble show slower 
damping than those from the majority DQDs but still 
they vanish on the time scales of several picoseconds. 



B. FWM response for long delays 

For longer delays, only F\ contributes to the signal. 
Out of the two components making up this term, F± + 
has the decay rate T + = f3 + T > T. This is the superra- 
diant component of the optical coherence, the decay rate 
of which reaches 2r when V ^S> A. The other, subradi- 
ant component iq_ has the decay rate T_ = f3-T < T, 
which vanishes in the limit of strongly coupled dots. The 
relative amplitude of these two components is f3 + //3_ , so 
that the superradiant one dominates if the coupling is 
strong, with a subradiant tail visible only at long delays. 

The decay of the FWM signal at long delay times is 
shown in Fig. [3Ja) for different strengths of coupling be- 
tween the dots. For V = 0, both decay rates are equal to 
r and one observes a usual exponential decay. As the cou- 
pling increases, the decay of the FWM response becomes 
non-exponential due to the presence of the sub- and su- 
perradiant components. For V ~ A, the decay sets off 
with an intermediate rate between T and 2r but it de- 
viates from an exponential form rather quickly. When 
V A the FWM response is almost completely dom- 
inated by the superradiant component, showing an ex- 
ponential decay with the rate 2r over a very long time 
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range, with only a weak tail corresponding to the subra- 
diant part of the signal. 

Since the decay rates r± depend on V/ft — V/(V 2 + 
A 2 ) 1 / 2 , which varies over the ensemble due to the in- 
homogeneity of A, they change from one DQD to an- 
other. In order to see how this inhomogeneity effect 
influences the nonlinear optical response we plot the 
time-integrated FWM response for a few values of cta 
in Fig. [3jb). The influence of this parameter is rather 
small and appears only for rather long times. The major 
effect is some softening of the shoulder which marks the 
transition from the dominating superradiant contribution 
to the subradiant tail. However, this happens only when 

cr A > n. 

Interestingly, as long as only the leading order response 
is considered, the biexciton shift appears only in the 
short-living term F3. At longer times, this kind of cou- 
pling does not affect the FWM response at all. Thus, 
there are no beats from molecular biexcitons in the FWM 
response. 



V. CONCLUSION 

Our results show that the shape and decay rate of the 
time-resolved FWM signal provide rich information on 
the properties of the DQDs in the ensemble, including 
coupling between the QDs. In the time- integrated sig- 
nal for long delay times, one observes a transition from 
the regime of independent decay (with the usual decay 
rate) to superradiant decay (double rate) via interme- 
diate cases of non-exponential decay. This transition is 
driven by the interplay of the energy mismatch between 
the QDs forming the DQD system and the coupling be- 
tween them and is only weakly affected by the distribu- 
tion of the energy mismatch. The superradiance effect 
observed in the decay of the FWM response is strongly 
stabilized by the coupling between the dots and can sur- 
vive even when the energy mismatch is many orders of 
magnitude larger than the radiative line width. 

Unlike ensembles of individual QDs, the DQD sam- 
ples show interesting and meaningful features also in the 
TI signal at short delays (several picoseconds). In this 
range of delay times, the signal has the form of oscilla- 
tions which result from optical beats between the radia- 
tion emitted from different DQDs. The form and decay 
time of these oscillations depends on the presence of a 
resonant subensemble (a subset of DQDs with matched 
transition energies) in the inhomogeneous ensemble. If 
such subensemble is present, its contribution dominates 
the short-delay response yielding a direct access to the 
strength of interaction between the dots forming the 
DQDs. 

In our study, we assumed that the surface density of 
DQDs in the ensemble is not very high and the distri- 
bution of the transition energies is rather broad so that 
collective effects on the level of the ensemble are absent. 
Including such effects would require applying a more gen- 



eral theor y 31 i 32 i 33 , as the DQDs are spaced by distances 
comparable or larger than the emitted wave length so 
that propagation and retardation effects would come into 
play. 
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APPENDIX A: THIRD ORDER RESPONSE: 
APPROXIMATIONS 

In this Appendix we give the full formulas for the com- 
ponents of the third order response and present the de- 
tails of the approximations that lead to Eqs. (I20al c). 
Substituting Eq. flT7|| into Eq. © one finds the FWM 
response signal in the form of Eqs. (| 18|) and (fT9|) . with 



*i±(A) = /^exp 



x exp 



(2iQ/h±(3 ± T) 2 T 2 c{ 
8(1-1 

-0±Tt 



P±T(t 



r) 



$2±(A) = exp 



2(l+r>|/^) 
(2ifi/fi± f3±T) 2 T? ei 



(Ala) 



8(1 

x exp [=F2if2/fir] 



a 2 E /h 2 ) 



'rcf u £ 



x exp 



-Tt 



P±T(to - t) 



2(1 



(Alb) 



$ 3 ±(A) = -/3 T exp 



{2iQ,/h±{2 + fe)T)\ 



2^2 i 

rel 



8(1 



x exp 



x exp 



x exp 



x exp 



Tl { v B (VB±2n) 
<4) 



ref"E 



2{H 2 

. (2 + /3 T )rr 2 f U 5 
'2(fi+r r >|A) 

~i(v B ±2n)T/h- 



iV B (t -T) 



-(l + /3 T )rr- 



(2 + f3 T )T(t -rj 



2(l + r>I/^) 



The values of the radiative decay rate T and the inho- 
mogeneous ensemble broadening of the transition ener- 
gies erg are of the order of [ieV and tens of meVs, respec- 
tively, while the typical duration of a reference pulse is 
about 100 fs or less. Based on these values, Eqs. (jAlal -c) 
can be considerably simplified. 

In the first exponents in Eqs. (|Ala( b). one can write 
(consistently with the approximation r <C Q/H used 
throughout this paper) 



{2in/h±p ± r) 2 T ] 



rcf 



m 2 T 2 /h z ±4i/3±rr 2 ef . 



Moreover, typically r-r rc f ~ 10 4 and flT re f/h ~ 0.1, 
hence the imaginary part can be safely neglected. The 
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same argument holds for the first exponential term in 
Eq. (jAlcp . so that in all three exponents this term can 
be replaced by exp[-(l/2)fi 2 /(^ 2 + T? e{ a%)]. 

Because of the Gaussian term in Eq. (flTj)) . it is clear 
that the measured signal is of considerable magnitude 
only when \to — t - | < h/uE, which reflects the 'photon 
echo' nature of the FWM response. Since hT/dE <C 1, 



the very last terms in Eqs. (|Alal -c). proportional to 
T(tQ— t), can be discarded. The third exponential term in 
Eq. (|Alc|) is also negligible, as r rc fr <C 1, while TrefVe/ft 
is typically of the order of 0.1 (for a biexciton shift of a 
few meV). With these approximations one arrives at the 
Eqs. pOak). 
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